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We present a detailed analysis of two-state trajectories obtained from force-clamp spec- 
troscopy (FCS) of reversibly bonded systems. FCS offers the unique possibility to vary 
the equilibrium constant in two-state kinetics, for instance the unfolding and refolding of 
biomolecules, over many orders of magnitude due to the force dependency of the respec- 
tive rates. We discuss two different kinds of counting statistics, the event- counting usually 
employed in the statistical analysis of two-state kinetics and additionally the so-called cycle- 
counting. While in the former case all transitions are counted, cycle- counting means that 
we focus on one type of transitions. This might be advantageous in particular if the equi- 
librium constant is much larger or much smaller than unity because in these situations the 
temporal resolution of the experimental setup might not allow to capture all transitions of 
an event- counting analysis. We discuss how an analysis of FCS data for complex systems 
exhibiting dynamic disorder might be performed yielding information about the detailed 
force- dependence of the transition rates and about the time scale of the dynamic disor- 
der. In addition, the question as to which extent the kinetic scheme can be viewed as a 
Markovian two-state model is discussed. 
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5} I. Introduction 

One of the major goals in biomolecular research is the understanding of how molecules such 
as nucleic acids and most prominently proteins fold into various conformations. Knowledge 

t— i of the underlying energy landscape is needed in order to understand how the large number 

>■ of possible conformations even a single chain can assume is funneled into a single stable 

S^ conformation. Subtleties of the energy landscape usually remained undiscovered and con- 

sequently folding is described in the framework of a simple two-state folder as the barrier 
crossing is the rate limiting step. 

With the advent of single-molecule experiments, however, the actual measurement of 
probability distributions came into reach that together with the tools of statistical me- 
chanics permitted to draw a more comprehensive picture of the folding and unfolding of 
biomolecules. For instance, single-molecule FRET (Forster resonance energy transfer) al- 
lows the separation of folded and unfolded subpopulations, which implies that folding and 
unfolding can be investigated under near native conditions pp. Besides optical methods, 
particularly mechanical single-molecule techniques have been successfully applied to study 
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the reversible folding and unfolding of molecules in an unprecedented way [2]. Here, addi- 
tion of a denaturant is not necessary since the molecules are stretched out of their native 
conformation along a predetermined coordinate. 

Even the simplest version of an energy landscape, two-states separated by a single en- 
ergy barrier, bears intricacies that typically escape ensemble measurements and sometimes 
even single-molecule approaches. The underlying multistate kinetics might be so involved 
that a single time series from, for example, photon counting statistics, is insufficient to un- 
ravel the complete mechanism. Single-molecule force spectroscopy in the reversible force- 
clamp mode, however, allows to deform the energy landscape by changing the externally 
applied force (3J 0|. Since the rate constants are strongly dependent on the external force, 
thermodynamics and kinetics can be tuned over a wide range, which essentially increases 
the information content that eventually might allow to access the multistate transition rates 
that are otherwise not amenable or at least provide a cue whether an unfolding/folding 
process is Markovian. 

Systems showing reversible bond-dynamics, such as certain biomolecules[51 El E] or 
specially designed molecules [8 J have been investigated experimentally by force spectroscopy 
in recent years and the impact of finite rebinding rates on the resulting force-spectra has 
been studied theoretically [9j [10l [11]. While in other fields of single molecule spectroscopy 
(SMS) [121 03] the analysis of trajectories has been used to gain information about the 
details of the underlying dynamical processes such an analysis has not been performed up 
to now in the field of FCS. 

The general theory of the treatment of the statistics of single molecule trajectories 
has been developed over the last decade with particular emphasis on reversible reaction 
schemes [TH [15] [16]. In particular, the analysis of the waiting time distributions (WTDs) 
can be used to distinguish between (non-Markovian) two-state behavior and more complex 
dynamics [T71 [18"] . In a previous paper [19], denoted as paper I in the following, two of us 
investigated the impact of non-Markovian dynamics on FCS experiments using a generating 
function approach introduced to the treatment of photon counting statistics by Brown [20]. 
We focussed on the mean number of transitions and its variance, the so-called Mandel 
parameter as a function of the external force. We showed, that similar to SMS, the Mandel 
parameter allows the detection of deviations from Markovian behavior. 

In the present paper, we apply the formalism developed in the quoted papers to trajec- 
tories as obtained from FCS experiments on systems exhibiting reversible bond-breaking 
dynamics. For example, the free energy landscape of the prominent GCN4 leucine zipper 
recently has been investigated in great detail as a function of external force using dual 
beam optical tweezers |21j. The idea of this paper is to provide a theoretical framework 
and feasible recepies that permit an in-depth analysis of (equilibrium) fluctuations as a 
function of externally applied force. In order to do so, we consider kinetic models with two 
'states' (ensembles) A and B (corresponding to 'bond-closed' and 'bond-open' states) and 
an arbitrary number of conformational substates within each of the A- and 5-ensembles. 
We will show that the possibility of changing the equilibrium constant for the transitions 
among the two ensembles is a key step to determine the detailed dependence of the phe- 



nomenological kinetic rates on the external force and thus to probe various models and 
particularly deviations from Markovian behavior. In favourable cases, the variation of 
the equilibrium constant might open the opportunity to determine the time scale of the 
dynamic fluctuations. 

The paper is organized as follows. In the next section we briefly review the standard the- 
oretical treatment of two-state trajectories and furthermore develop the techniques needed 
for our description of FCS. In particular, we compare the information that can be extracted 
from different counting protocols. We then propose a possible route of analyzing experi- 
mental (or simulated) trajectories and show how relevant information about the kinetics 
in the system investigated can be obtained. The corresponding calculations are performed 
for a simple model with two conformations in each ensemble (two-configuration model or 
two-channel model, TCM). 

II. Theory of event-counting statistics with application 
to FCS 

The general formalism for the analysis of two-state trajectories as observed in single- 
molecule experiments has been developed mainly in the past decade [HI [TS1 [22]- In this 
Section, we briefly review those formal techniques that we will use to analyze the trajecto- 
ries occuring in FCS of reversible systems and additionally discuss the theoretical treatment 
of different techniques of counting. 

The simplest (coarse-grained) model of reversible bond-breaking is a two-state model 
for the populations of the 'states' A and B, defined by the minimum regions in an energy 
landscape as sketched in FigjIJ If described as a Markov process, the master equation [23] 
reads: 

d t px(t) = -k x px{t) + k Y p Y {t) with X, Y = A, B and Y + X (1) 

where Isa = k(A — > B) and ks = k(B — > A) are the corresponding transition rates. In 
FCS, a constant external force f ex t is applied to the system under consideration (a single 
molecule or an adhesion cluster, for example). The dependence of the transition rates in 
an energy landscape with two (or more) minima on the external force can be approximated 
by Kramers theory, cf. for instance the treatment in ref.pj]. Quite generally, the applied 
force changes the apparent activation energy and therefore one has a strong dependence 
of the kinetic rates on the external force, cf. FigjTJ The activation energy for the A — >■ B- 
transitions will decrease as a function of the applied force f ex t whereas the one for the reverse 
transition will increase, meaning that £u(f e xt) increases and fcs(f ex t) decreases. Therefore, 
the equilibrium constant 

#(fext) = MfextVMfext) (2) 

strongly decreases as a function of f ext . As already noted in the Introduction, this fact 
allows to investigate vastly different physical situations by FCS that otherwise are only 
accessible via e.g. temperature variations in favourable cases. In the following, we will 
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Figure 1: Sketch of the free energy G(q, f ext ) as a function of the reaction coordinate defined 
as the pulling direction for various values of the applied force. For larger applied force, the 
right 'B' minimum becomes deeper (dashed and dotted lines) and the equilibrium constant 
A'(fe X t) decreases. 

always use the shorthand notation K = K(i ext ) and have in mind that K can easily be 
varied via variation of f ex t- Similarly, all rates are to be understood as depending on the 
external force. From the stationary solution of eq.pl), 
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it is evident that the system predominantly resides in one of the states if K strongly 
deviates from unity and that the reversible nature is most pronounced for K ~ 1. 



A. Statistics of transitions - a brief review 

If the statistical nature of the transitions between the A- and -B-states is important as is the 
case when dealing with trajectories, it is not sufficient to consider the global populations 
Px(t) as it is done in the master equation, eq.Q. Instead, one has to treat the individual 
transitions and describe their statistical properties. There are different ways to do so. 
One possibility is to consider a two-component Markov process consisting of the 'reaction 
coordinate' and the 'number of transitions' [T9| 120] . Another way of treating the statistics 
is to decompose the system into two 'half-reactions' describing the decay of either A or B 
and the full kinetic scheme is then obtained by 'adding' the half-reaction schemes p3]. The 
most general procedure consists in rendering the observed transition irreversible [15J. The 
two ways just mentioned can be viewed as special cases of this general method. 

In the following, we assume that we do not only have a single state for A and B 
but a multitude of (conformational) substates defining the A-ensemble and similar for 
the .B-ensemble. The populations of these ensembles are then described by the vectors 
P^ = (pa,i,Pa,2, ■ ■ ■ ,Pa,n a ) and similar for P^. The dimensions Nx give the number of 
substates in the X-ensemble, X = A, B. Of course, continuous models can be treated in a 
similar way. Our resulting kinetic scheme then reads as: 
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Here, the elements of the matrices Kj are the inter-ensemble transition rates ky,a;X,a' for 
each (X, a') — > (Y, a) transition: 

(K x ) aa , = k Y , a . x , a > [X ± Y) 

(K' x ) a>a > = S a , a 'k X;a with k x . a = ^J k Y , a ';X,a (5) 

a' 

Furthermore, the IIx consist of the rates for the transitions (X, a') — > (X, a) within the 
X-ensembles, the so-called exchange rates , ~f X - a ,a' characterizing the dynamic disorder: 

( n *)a,a' = -$a,a'JX;a + (1 ~ ^a, a ')lX;a,a' With 7x ;a = ^ 7X;«',a (6) 

Of course, the fact that the II x are master operators and therefore the (conformational) 
dynamics within a given ensemble is Markovian can be relaxed. The corresponding Green's 
function obeys: 



II , - K' A Ki 

-a n s — xv B 



G(t) = wc(t) with w=( a Ka A n K' ) (7) 



with the initial condition G(0) = E, where E denotes the unit matrix. 

In order to extract the observed transitions from the remaining unobservable ones, one 
decomposes the transition-rate matrix according to|15j: 

W = W + V (8) 

Here, W consists of the rates for unobserved transitions and the elements of V are the 
rates for the observed transitions, cf. the discussion below. Note that the system described 
by the transition-rate matrix W is irreversible, i.e. there is no finite stationary population 
of the states. For the Green's function of this irreversible system one has: 

G'(£) = W'G'(t) or G'(t) = exp (W't) (9) 

As already indicated above, in the present paper, we will consider two different choices 
for the observable transitions, i.e. for the matrix decomposition in eq.pj), cf. FigJ2J 

One choice is the one that was used in our former treatment in I and consists in the 
consideration of whole cycles: 

U A -K A \ ,, _ / K B 



Here, only the (B — > v4)-transitions are counted as observable, as shown in the left part of 
FigfJJ Of course, all results obtained from this choice are identical to the corresponding 
ones when interchanging A and B. In the following, we will denote this type of observing 
only full cycles as 'cycle-counting'. 

The other type of decomposition is the one that typically is employed when the statis- 
tical properties of the WTDs (to be defined below) are considered. It consists in counting 
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Figure 2: Example trajectories schematically showing the two different ways of counting 
as indicated by the arrows. Left: cycle-counting (only B — > A-transitions are counted); 
Right: event-counting (all transitions are counted) 



all inter-ensemble transitions and will be called 'event-counting' in what follows, cf. the 
right part of Fig|2] 
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Here, the underlying system is that of two half-reactions [13] and one has from eq.Q 
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G' B (t) 
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(12) 



For this decomposition several indicators of complex dynamics have been introduced and 
discussed in the context of SMS [21 ESI \M EZ] • 

An important quantity that can be extracted from single molecule trajectories and that 
also will be of interest for us is the distribution of the time between consecutive observable 
transitions [To]. 

V(t) = l T VG / (r)VP st '/(n) (13) 

Here, 1 T = (1, ■ ■ ■ , 1) is a summation row vector and the stationary populations obey 
WP st ' = and are normalized, l T P st = 1. Furthermore, we defined the mean number of 
observable transitions per unit time: 



St. 



(14) 



(n) = 1 T VP 

Here and throughout the paper we use the stationary flux VP st, /(n) as initial condition[T 
because this can easily be realized by starting the counting from any event of a represen 
tative trajectory [T41 IT5] . 



In case of cycle-counting the distribution of transition times, eq.(13), corresponds to 
the 'turn-over time distribution': 



V c (t) = l T V c G / c (t)V c P st -/(^)c with (n) c = l T V c P st (15) 

On the other hand, for event-counting, one usually considers the WTDs[H], given by: 



$ x (t) = l r Y K x G' x (t)K Y Pf/(n) x with (n) 



x 



l^KyPp- 



Y^X = A,B (16) 



Here, we defined the summation row vector 1^ with dimension N X - In order to quantify 
correlations in the trajectory, one studies the distribution of consecutive waiting times: 



$x,y(ti,£ 2 ) = l T x K Y G Y (t 2 )K x G' x {t 1 )K Y Pp-/(n) 
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If the events are uncorrelated, one has $x,y(^i ? ^2) = &x{ti)& Y {t2) for all combinations of 
X and Y and the kinetic scheme is called reducible, meaning that it can be reduced to a 
two-state model, albeit usually with non-exponential WTDs[18j. 

Another quantity of interest is the probability to observe a given number of transition 
in a prescribed time interval. The generating function for the moments of this distribution 
can be written in a compact notation as[T5l 



F(z,£) = l T e (w ' +zV)f VP s 7(n) (18) 

and the moments are computed in the standard way from F(z,t) via differentiation, for 
instance: 
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and (N(N-l)(t)) 



d 2 F(z } t) 



dz 2 



from which one can calculate the Mandel parameter 
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(19) 



(20) 



It can be shown that both, (N(t)) and Q(t), are related to certain event-correlation 
functions [15J, but we will not dwell on this further in the present context. We will use 



the expressions (19) and (20) for both, event-counting and cycle-counting in the following 



section. This means that we use the respective decomposition, eq.(10) or eq.(ll), in the 



general expression for the generating function, eq.(18) 



B. Event-counting versus cycle-counting 



As is to be expected from their definition, eq.(15) and eq.(17), there exists a relation 



between the turn-over time distribution V c (t) and the distribution of consecutive waiting 
times. It is shown in Appendix A that this relation is given by: 
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where we used that (n) c = 1 T V C P S< - = I^KbP^ = (ti)a- The system starts in state 
B before a transition to A takes place. In A, the system evolves for some time until it 
goes over to B and the final B — y A-transition is counted. In the cycle-counting scheme 



all intermediate transition-times are summed yielding the integral in eq.(21). If the A —y 



S-transitions are counted as observable, one has to interchange A and B in the above 
expression. 

As moments of the WTDs may be more easily extracted from experimental (or simula- 
tion) data than the WTDs themselves |1 8], we give the expressions for these moments. The 



derivation of these expressions is presented in Appendix A. There the general expressions 



for the nth moments are given in eq.(|A.ll) from which one obtains: 

(r 2 ) x = 2.i x G' x (0)Pf/(n) x 
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Here, G^-(O) denotes lim^o G' x (s) with G' x (s) = f^°dtG' x (t)e st being the Laplace trans- 
form of G' x (t). 

We mention that the form of the exchange matrix Tl x given in eq.([6]) yields ljllx = 
and from this one can show that i^KyPff = ly] 



llK x Pf and thus 



(n) A = (n) 



B 



(23) 



Only in generic non-equilibrium systems we expect this relation to be violated. The first 
moments of the WTDs can be used to define phenomenological rate equations. We define 
the effective transition rates [25] 

k x = l/(r) x (24) 



and find, using the expression for (r) x , eq.(22), and eq.(23): 



k A pf = k B pf with Px - = 1 T X P X - 



(25) 



It is thus meaningful to propose the rate equation given in eq. (fTl) for the dynamical evolu- 
tion of the system on a coarse-grained level, cf. the following Section. 



For the moments of $x,y(^i,^2) one finds from eq.(A.13) 
(r) XY = l r Y G' Y (0)K x G' x (0)Pf/(n) 
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Reducible schemes can be identified via the factorization property 

(r n )xv = (r n )x(r n )Y 



(26) 



(27) 



Furthermore, we can give the moments of the turn-over time distribution V c (t). The 
actual calculation is presented in Appendix A and here we give the results: 



(t) c = (t) a +(t) b 
(t 2 ) c = (t 2 ) a + (t 2 } b + 2(t} 



AB 



(28) 



The first moment in case of cycle-counting equals the sum of the moments of the WTDs 
relevant for event-counting and correlations between the transitions only affect the second 
(and higher) moments. In particular, cycle-counting can be used in order to investigate 
the question whether the system under consideration can be described by reducible kinetic 
scheme in favourable cases. In order to decide this, one considers the variance, (St 2 ) c = 
(^ 2 )c - (t) 2 c which is given by (5t 2 } c = (5t 2 ) a + (5t 2 ) b + 2((r} AB - (t) a (t)b). For a 
reducible scheme the last term vanishes indicating that the transitions are uncorrelated. 
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III. Possible analysis of trajectories 

In order to apply the theoretical results presented one has to rely on model calculations 
in most cases because generally valid results are rare. In the following we will discuss 
a possible way for the analysis of two-state trajectories as they can be obtained from 
FCS-experiments on reversibly bonded systems. As mentioned above, one distinguishing 
feature of FCS is provided by the possibility of varying the equilibrium constant K by 
a large amount allowing us to discuss details of the underlying kinetic scheme that are 
difficult to access by other means. The particularly interesting question of deviations from 
Markovian/Poissonian behavior of the system under study will be treated via considering 
the impact of dynamic disorder [28J in terms of so-called exchange models. 

Before we proceed, let us briefly discuss possible origins for dynamic disorder. One 
obvious reason for the existence of a large number of transition rates is given by the fact 
that in many biomolecules an ensemble (e.g. folded or unfolded) consists of a huge number 
of possible molecular conformations. Therefore, conformational changes within a given 
ensemble play the role of exchange processes taking place in addition to the A — > B- 
transitions. Another possibility that has already been discussed in the literature is given 
by fluctuations in the transition state ensemble [211 ED] and this has also been observed 
experimentally pTj . In terms of the simplest conceivable model, the Bell model, one has a 
dependence of the form k Ya -Xa' ~ e _axfcxt , where ax denotes the distance to the transition 



state, see eq.(29) below [52^. A distribution of the latter thus results in a corresponding 
distribution of k x ,a t Y,a'- Sofar, the concept of these 'bond- heterogeneities' has been ap- 
plied to systems showing irreversible rupture events. An analysis of reversibly bonded (or 
refolding) systems opens the opportunity to study the impact of such fluctuations in the 
transition state ensemble in much more detail. If the ax-fluctuations take place on the 
experimental time scale, the system exhibits the features of dynamic disorder. 

In the present paper, we will mainly use the TCM with only two substates in each 
ensemble for illustrative calculations but we try to keep the discussion as general as pos- 
sible. Details of the calculations concerning the TCM are presented in Appendix C for 
convenience of the reader. Before we start considering the impact of dynamic disorder, we 
propose to extract averaged information from the trajectories. 



A. Determination of force-dependent mean kinetic rates 

We suggest that the first step in the analysis of two-state trajectories consists in the deter- 
mination of the first moments of the WTDs. This should be possible with high statistical 
significance and allows to define the effective transitions rates kx on a phenomenological 
level, kx = (t)x 1 °f- eq.(24)[25j. (We will simply write k x instead of kx in the remainder 
of this section.) The resulting phenomenological two-state model consists of only a single 
A- and a single 5-state and the temporal evolution of the populations of these states is 



governed by the master equation, eq.Q. 



As pointed out several times, a broad range of kinetic rates fcx(fext) is accessible via 
FCS. This is because already in the phenomenological Bell model[32] one has an exponential 
dependence of the kinetic rates on the external force, 



M f ext) = k A (0)e 



aiis: 



and fee (f e 
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(29) 



where ax are the distances from the transition state location scaled by the inverse temper- 
ature /3 = 1/T (we set the Boltzmann constant to unity). There are thus two parameters 
that are to be determined from the (t)x as a function of the external force. Note that in 
the Bell model one neglects any f ext -dependence of the positions of the free energy minima 
and the transition state and also higher order powers of f ext [5j [191 E21 EI] . In the present 



paper, we will use the simple form given in eq.(29) for model calculations and all times 
will be given in units of fc^(O), cf. FigJ3^a). Because all of our arguments also apply to 
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Figure 3: (a): kinetic rates k A = (t)^ 1 and ks = (t)b 1 as a function of the external 
force (top) or of the equilibrium constant K(i cxt ) (bottom). The parameters are chosen 
as a a/ 'ft = 0.3nm and a^jft = 0.7nm and the temperature is ft^ 1 = 300-K" = A.lApNnm. 
Furthermore, k A (0) = Is -1 and fcs(0) = 10 4 /c^(0). (b): Typical trajectories of the two- 
state model for K — 1 and K — 10 obtained from kinetic Monte Carlo simulations for the 
same parameters. 

other single-molecule techniques and are not restricted to an analysis of FCS-data, in the 
following we will always refer to the equilibrium constant instead of the external force. Of 



course, the simple form for the rates given in eq.(29) can be relaxed but at the expense of 
introducing more parameters. 

Experimentally, it should be possible to determine the f ex t-dependence of the kinetic 
rates for a rather large range of forces. In particular, we propose that it is possible to 
quantify deviations from the simple Bell model without relying on some specific model 
for the f ex t-dependence. This fact distinguishes FCS from standard applications of force 
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spectroscopy using a linear force ramp [35] . where a model- free analysis appears to be 
difficult [33]. 

In this context, the importance of the opportunity to apply different counting schemes 
becomes apparent. When considering trajectories with different values of the equilibrium 
constant, as shown in Fig|3[b) for K = 1 and K = 10, it is obvious that for a finite 
temporal resolution an analysis in terms of event-counting becomes difficult for situations 
in which K differs strongly from unity. In such cases, cycle-counting appears better suited 
and can also be used to determine (r) x as for K ^> 1 one has (r) c = (r) A and for K <C 1 



(r) c = (t)b holds, cf. eq.(28). Furthermore, for values of K that allow an analysis in terms 
of both counting schemes a check for internal consistency of the data analysis is given by 
the general relation (r) c = (t)a + ( t )b- 

An independent test of the analysis performed can be obtained via a comparison of the 
values for kx determined from the first moments of the WTDs with the mean number of 
transitions in a given time t, (N(t)). The calculation of the full time- dependent expressions 



for (N(t))^ ev ^ and (N(t))^ are outlined in Appendix B and the results are given in eq.(B.3). 
From these expressions one can see that in both cases after a time on the order of the inverse 
relaxation rate (k A + &b) a behavior linear in time is reached: 

(N(t))^ = 2(N(t))^ = (n)t ; t » (k A + ks)' 1 (30) 



Here, the mean number of observable transitions per unit time according to eq.(14) is 

((n)^ = (n) A + (n) B ): 

<»> s <»>'-' = 2 <"> w = 0t> (31) 

As is to be expected, in case of event-counting one has twice the number of transitions as 
compared to cycle-counting. 

From the above analysis it should be possible to obtain the detailed dependence of 
the phenomenological kinetic rates on the external force. Therefore, one can analyse this 
force-dependence in more detail and for instance use existing models in order to extract 
information about the location of transition states in the free energy landscape and details 
of their properties. 

As a second step in the analysis of two-state trajectories we propose to check for de- 
viations from simple Markovian two-state behavior of the kinetic scheme under study. 
As already discussed in I, the Mandel parameter is well suited to detect deviations from 
Markovian behavior. The Q-parameters Q(t)^ ev ^ and Q(tp c ^ are time-dependent on a scale 
of (kA + ks) and then approach their constant long-time limits Qoo and Qoo . The calcu- 
lation closely follows the corresponding one of (N(t)) and one finds, using the definition 



eq.(20): 



In Figj4j we show Q^ and Qoo as a function of the equilibrium constant and thus of the 
external force, cf. Fig|3] By definition of the Mandel parameter, one has deviations from 
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Figure 4: QS? (red, upper line) and Q£J (black, lower line) as a function of the equilibrium 
constant. The parameters are the same as in FigJHj The full lines are the analytical 
expressions, eq.(32), and the points are results of a kinetic Monte Carlo simulation. 



Poissonian behavior whenever Q ^ 0. It is shown in Appendix B that in case of event- 
counting one has Poissonian behavior for K = 1 and for cycle-counting one finds a Poisson 
distribution for K <C 1 or K ^> 1 . Because for the determination of the moments of the 
event number probability one has to count all relevant transitions, we expect that also 
for the determination of the Q-parameters cycle-counting will become the more reliable 
method for K strongly deviating from unity. The information content, however, is the 
same for both counting schemes. 

Non-Markovian effects will of course also be reflected in the WTDs. In particular, 
deviations from exponentiality indicate complex kinetic behavior. As has been noted 
by Flomenbom et al.[T8], experimentally it usually will be easier to determine the mo- 
ments of the WTDs quite accurate whereas the the WTDs will often be very noisy. As 
shown in Appendix B, for the phenomenological two-state model, the WTDs are exponen- 
tial distributions, $>x(t) = kxe~ kxt , and the corresponding moments are trivially related, 
(t 2 )x = 2(t)x- Furthermore, the kinetic scheme is reducible, $x,y(^i,^2) = < ^ > x(^i)-$y(^2)- 



If the Q-parameter as a function of the equilibrium constant follows eq.(32) and the 
second moment of the WTDs obey (t 2 )x = 2(r)^-, one can safely conclude that the system 
under study can be described by a simple kinetic two-state model and the analysis is 
complete. 

B. Analysis of complex kinetic schemes 

If the above analysis of the two-state trajectories indicates deviations from Markovian 
behavior no general procedure of quantifying these can be proposed. This is because the 
analysis of two-state trajectories does not allow to determine the underlying kinetic scheme 
unambiguously [16j [TTJ, [18]. The impact of dynamic disorder or dynamic heterogeneities 
can thus only be discussed in terms of specific models and different indicators aiming at 
quantifying the relevant time scale have been introduced J2U ESI |26l [27] . When treating the 
effect of conformational changes or dynamic disorder via the introduction of exchange rates 
lx-,a,a' m eq.([6]), one can distinguish two limiting dynamical regimes. For exchange rates 
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large compared to the inter-ensemble transition rates, the so-called motional narrowing 
regime, one obtains the phenomeno logical two-state model just discussed. In the other 
limiting case of vanishing 7x ; a.o' the relaxation properties are determined completely by 
the distribution of the kx, a -Y,a' (static disorder regime). In both cases, the time scale of 
the dynamic disorder cannot be determined from the data. 

If deviations from Markovian two-state behavior are observed, there are two questions 
to be posed. First, one would like to know the time scale of the dynamic disorder. Fur- 
thermore, it is interesting to determine whether the system still is reducible meaning that 
one has two-state behavior. This latter point is of particular interest when dealing with 
biomolecules like proteins, where two-state behavior is not always quantified easily. 

Time scale of fluctuations 

One limiting scenario that should clearly be observable in the data is the case in which the 
time-dependent Q-parameter does not reach a constant limit Qoo at long times. As this 
limit is reached on the time scale of the slowest rate in the system, one can conclude that 
on the experimental time scale one is in the static disorder limit. A further analysis of 
the data for the Mandel parameter usually will be difficult because numerical calculations 



show that the linear ^-dependence of Q(t) at long times as given in eq.(C8) for the TCM 
will only be reached on a time scale much longer than the inverse inter-ensemble rates. 

In this limit of static disorder, one can give general expressions for the moments of the 
WTDs. Defining the average of quantities in the A-ensemble via 

N x 

{A) x = Y, A «Pt« ( 33 ) 



a=l 



one finds for instance 



Nx 

(n) X = l£KyP£ = Y, k X;aPx]a = (k X ) X 
a=l 



where kx ;a denotes the overall escape rate from the substate (X,a), cf. eq.(J5l). The first 
moment of the WTDs follows from eq. ( 24 ) : 



r)x = k x 1 =pf/{k x )x (34) 



with p s x = i^Pjf, cf. eq.(25). For the second moment, (r )x, one can give an analytical 



expression because in the static disorder limit one has [G^-(0)] aiK = S aiK k x . a . Using this 



expression, one finds from eq.(22) that 

(r 2 )^ = 2(kx) x 1 (k?) x (35) 

The analysis of the mean rates discussed above yields information about (kx)x an d one 
can try to obtain the width of the underlying distribution of rates from (k x )x- For 
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instance, for the TCM, one has {k x )x — P x '(k x ) with p s x given in eq.pl) and one finds 
(k x 1 ) = (kx)/[(k x ) 2 -(5k x )]. 

Another example that can be used to illustrate the relation between the differ- 
ent moments is the model introduced by Raible et al. mentioned abovej29j [30] . Us- 
ing the Bell expression given in eq.(29), kA(aA,iext) = kA{0)e aAicxt and a Gaussian 
distribution of the form g(au) = N' 1 exp [— (a^ — (a J 4)) 2 /(2A^)]9(a J 4) one finds that 
(k^x = Px' J^°daAg(o:A)kA(o:A,iext) n - The deviations from Markovian behavior are best 



characterized in terms of the ratio (Sr 2 ) x / ( T )x = ( r2 )x/( r )x ~ 1 which equals unity in 
the Markovian limit. In the static disorder limit relevant here one finds 



(5T 2 )x/(T} x = 2(kx)x(k x 1 ) 



x 



In Fig|5] we show this quantity as a function of the width of the a^-distribution. It is 




A A / <a A > 

Figure 5: (5t 2 )a/(t)\ for a model of bond-heterogeneities in the static disorder limit [29| 
30] as a function of the width A^ of the Gaussian distribution of the distance to the 



transition state, g{a.A) 
0.3nm. 



N- 



-i 



exp [— {a a — (a y i)) 2 /(2A^)]0(a J 4). Here, we used {a a) / ft 



obvious, that independent of the equilibrium constant the Markovian limit is reached for 
small widths A a- For larger widths a change in the equilibrium constant could help to give 
an estimate for the width of the distribution. (When applied to experimental data, values 
Aa ~ (0.5 ■ ■ ■ I) (a a) have been found [30].) 

A possible cross-over from the static disorder limit to dynamic disorder on the time 
scale of the experiment /simulation is most easily observed when considering Q(t). If the 
time scale of the dynamic disorder is reached in the time window used in the analysis, Q(t) 
is known to reach a finite ^-independent value at long times. This limit Qoo is shown in 
Fig J6] for the TCM for various values of the exchange rates ■ja and 75. It appears that it is 
mainly the smaller exchange rate that determines the overall behavior and also the obvious 
deviations from the Markovian value (dotted lines). As we have already noted in I, the 
determination of the Mandel parameter allows the 'detection' of deviations from Markovian 
behavior but this does not allow to discriminate among different kinetic schemes. However, 
a variation of the equilibrium constant might be used to attempt a determination of the 
mean exchange rate. 
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Figure 6: Q^ versus external force for various exchange rates ja and 7# for the TCM 
defined in Appendix C. The parameters chosen are fc^, a (f ex t) = 1.99/c(f cxt ), &4,b(fc X t) = 
0.01fc(f cxt ), where fc(f ex t) = kA(0)e aAicxt as defined in eq.(29). This means that according 



to eq.(C3) one has (kx) = k and ax = 0.995/c. The other parameters are those given in 



the caption to Figj3 The dotted red line is the Markovian limit given in eq.(32). 



We have pointed out above that deviations from Markovian two-state behavior also 
manifest themselves in the deviations of the moments of the WTDs from the simple be- 
havior (t 2 ) x = 2(t) x - I n FigJ7]we show {St 2 ) a /(t} 2 a for the TCM as a function of the 



exchange rate. (The expressions for the relevant moments are given in eq.(C.ll).) In the 




Y A / < k A> 

Figure 7: (St 2 )a/(t) a versus ja/ ' {k A ) for the TCM showing the transition from the static 
disorder limit for small 7a to the motional narrowing regime for large 7a- The upper right 
inset shows the dependence on the 'width of the rate distribution', a a- Here, the values used 
are 0.995: (A; A>a (f cxt ) = 1.99fc(f ext ), k A , b (f ext ) = 0.0l£(f ext )); 0.949: (_^, (f ex t) = l-9fc(f ext ), 
&A,b(fext) = 0.1fc(f ext )); 0.707: (Au,a(fext) = l-5/c(f ext ), k A ,b(f ex t) = 0.5A;(f ext )). The left inset 
shows (St 2 ) a/ (t) a as a function of the equilibrium constant indicating that in some cases 
the cross-over between different regimes might be observable. The remaining parameters 
are the same as in Figj6j 

static disorder limit one has the expression given above and in the limit of motional nar- 
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rowing, the model reduces to the Markovian two-state model and one has (St 2 )a = (t)\. 
The dependence on the width parameter is shown in the inset. Also shown is the depen- 
dence on the external force in terms of the equilibrium constant. From this plot it becomes 
apparent that in favourable cases one should be able to observe a transition from static to 
dynamic disorder or from dynamic disorder to motional narrowing. For the specific model 
considered here, we find that one has about three to four orders of magnitude of exchange 
rates which define the regime of dynamic disorder. In this regime, one might attempt to 
determine the exchange rates using one of the methods that have been discussed by Cao 
and coworkers P H UE \ [251 [27] . Additionally, the long-time limit of the Q-parameter can be 
helpful to determine the mean time scale of the fluctuations. 

Reducibility of the kinetic scheme 

We have mentioned already that the question whether the kinetic scheme under study 
can be considered as a two-state system is of interest in many applications, in particular 
when concerned with the folding dynamics of biomolecules. As discussed in the context 



of eq.(17), a kinetic scheme is reducible, if the joint probability function &x,Y(ti,t 2 ) of 
successive waiting times factorizes, $x,y(^i,^2) = ^xit^^yih), for all combinations X, 
Y. We note that $x,y(^i 5 ^2) contains all information about correlations in the trajectory 
and in principle is obtained from the data by constructing a histogram of the intersections 
of successive (X — y Y~)-transitions followed by (Y — y A)-transitions[16j. In ref.[18] the 
various connectivities giving rise to a reducible scheme are discussed. In Appendix A, we 
show that in the static disorder limit also models with a global connectivity among the 
substates of the two ensembles can be reducible for a certain form of the transition rates, 



cf. eq.(A.17) 



In order to test experimental data for reducibility, in general one has to perform a 
detailed comparison of the WTDs $a(£), $b(£) and $x,y(^i,^2) or the corresponding mo- 
ments. In particular, the determination of $x,y(^i,^2) (or the corresponding moments 
(T n )xy) from event-counting statistics might be difficult with the required accuracy. In 
this case, cycle-counting can provide the desired information with much higher significance. 
A comparison of the variances (St 2 ) c and {5t 2 )a, (5t 2 )b yields important information about 
( t )ab, cf. the discussion at the end of Section II. 

We show the behavior of (5t)ab = {t)ab — (t) a{t)b for the TCM as a function of the 
exchange rate in FigM cf. eq.(|C.14). It is evident that the system becomes reducible for 



large exchange rates, i.e. the motional narrowing regime. The transition to this reducible 
regime depends on the ratio of the exchange rates in the two ensembles (upper panel) and 
also on the equilibrium constant (lower panel). In some favourable situations it might be 
possible to observe the transition to the reducible regime via variation of the external force. 
It is obvious from the expression for (St)ab that the maximum value increases with larger 
ax- Apart from the motional narrowing limit, the TCM becomes reducible in the limit 
^A,b — )• in which case the kinetic scheme consists of a single 'gateway state' [TS] . 

Another situation resulting in a reducible scheme that might be of experimental rele- 
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Figure 8: (St) ab / ((t) a(t) b) versus 7a/(&a) for the TCM. Upper panel: different exchange 
rates are assumed for the A- and the 5-ensemble, ^b/ia = 0.01, 0.1, 1, 10, 100 from right 
to left. The transition to the reducible regime depends on their ratio. The red line is for 
1b = 1a- Lower panel: transition to reducibility as a function of the equilibrium constant. 
The parameters are the same as in Fig|6j 



vance is given by the limit of one exchange rate ^x being much larger than all other rates, 
cf. FigJ8J (upper panel). This scenario might apply if the X-ensemble represents a stiff 
bonded system (with many closed hydrogen bonds) and in the other ensemble the bonds 
are opened and therefore there is more structural flexibility for conformational fluctuations. 
The dynamics of the escape from the X-ensemble is a Poisson process but the dynamic 
disorder present in Y(=£ X) is reflected in the measurable quantities like Qc» (cf. Fig, 6). 



V. Conclusions 

We have investigated the statistics of two-state trajectories as they are expected to occur 
in force-clamp spectroscopy of reversibly bonded systems. As has been shown earlier, one 
can detect the deviations from Markovian behavior from an analysis of the moments of the 
distribution of the overall number of events. In particular, the Mandel parameter allows the 
direct observation of dynamic disorder. We have discussed two different ways of analyzing 
the trajectories. One can either count all transitions, denoted as event-counting, or one 
counts only one kind of transitions ((A — > B) or (B — > A)), termed cycle-counting. We 
have shown the interrelation between the WTDs and their first two moments for the two 
counting schemes. Depending on the value of the equilibrium constant one of the methods 
might be superior to the other if applied to real data with finite temporal resolution. In 
particular, for values of the equilibrium constant strongly deviating from unity, the system 
mainly resides in one of the states (either A or B) and event-counting might not be feasible. 
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In this situation, cycle-counting still gives access to most dynamical information required 
to analyse the data. 

We propose that the first moments of the WTDs as a function of the external force 
can be used for an analysis of the dependence of the kinetic rates on the external force. In 
particular, one does not rely on fitting explicit expressions derived from some model to the 
data. The consisteny of the analysis can be checked by comparing the phenomenological 
transition rates obtained from the first moments of the WTDs with the mean number of 
events. 

An analysis beyond the determination of the mean rates can be proposed and allows 
to determine whether the kinetic scheme is determined by dynamic disorder or if one is in 
one of the limiting situations of static disorder or motional narrowing. In order to discuss 
these issues we have performed model calculations for one particular kinetic scheme, a two- 
configuration (two-channel) model, but we expect the general features not to depend too 
sensitively on the model used. In particular, the Mandel parameter should behave quite 
different in the mentioned limits. Thus, its determination gives a first hint on the time 
scale of dynamic fluctuations. 

As the various moments depend on the inter-ensembles transition rates and the ex- 
change rates in different ways, we expect that a comparison of the first and second moments 
of the WTDs as a function of the force dependent equilibrium constant can be helpful in the 
determination of the average exchange rate. This holds in particular, if the indicators of dy- 
namic disorder that have been discussed in the past are applied additionally [!Ml |25| |26 | [27] . 
Also the long-time limit of the Mandel parameter can be used for the purpose of deter- 
mining the fluctuation time scale. It is a major advantage of FCS to provide a means of 
varying the equilibrium constant by a huge amount due to its dependence on the external 
force. Therefore, FCS might be superior to other methods in determining the life-time of 
the dynamic disorder. 

Even though a distinction of different kinetic schemes is not possible in general, con- 
sideration of the moments of the WTDs for event-counting and cycle-counting might in 
favourable cases allow to determine whether the kinetic scheme underlying the observed 
two-state dynamics is reducible or not. This analysis might be very helpful in the discussion 
of the folding pathways in biomolecules. 

In summary, we propose to analyze two-state trajectories as they can be determined 
from FCS experiments in order to gain information about the force-dependence of the 
kinetic rates, the deviations from Markovian behavior and about the time scale of dynamic 
disorder, if present. We hope that our analysis will be helpful in the design of explorative 
experiments in the future. 
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Appendix A: Some properties of the WTDs 

In this Appendix, we derive the relation between the WTDs and the turn-over times 
distribution and give the general expressions for the moments of the WTDs. 

Relation between V c (t) and Qx,Y(ti,h) 



We start from the Laplace transform of eq.(15), 

V c (s) = l T V c G' c (s)V c P s 7(n) c 



(A.1) 



and write, using the definition for W^., eq. (JlO]) 

G'M = 
where we used the abbreviation: 



sE A - W^ 

-K A sE B - W' B 



W, 



X 



n 



X 



K' 



x 



(A.2) 



Recognizing that G' x (s) 1 
block matrix, one finds: 



[s~E x — W^-) and using the formula for the inversion of a 



G'Js) 



G' A (s) 

G' B (s)K A G' A (s) G' B (s) 



Inserting this decomposition into eq.(A.l) yields 

V c (s) = 1 t a K b G' b (8)K a G' a (s)K b I>£/{ti)a 



(A.3) 



because of (n) c = l T V c P st ' = (n) A . If after inverse Laplace transform this expression is 



compared to the definition of $a,_b(^i,^2) in eq.(17) one obtains eq.(21). 



Moments of the WTDS 

The moments of <&x{t) are defined in the following way: 

(r n )x= dtt n $ x (t) and (r n ) XY = / dh / dt 2 «^,y(ti,t 2 ) 

Jo Jo Jo 



We start from the definition of the WTD, eq.(16), and use the following identities 



1 X U X = and l^KyAy = l^K y A y 



(A.4) 



(A.5) 



where Ay is an arbitrary vector of dimension Ny- We will furthermore use the detailed 
balance condition 

WP st = -> W x Pf = -KyPf* (A.6) 



19 



where P s *' = (P^', P s b) t and the second identity follows from the structure of the transition 



rate matrix W' ei , given in eq.(ll) and additionally we used the definition of W^-, eq.(A.2). 
We now use these expressions to slightly rewrite eq.(16) for & x (t). In the expression 
l Y K x G' x (t)K Y Pf, we use eq. gjHb in the form 1^K X A X = 1^K^A X = 
-i^W^A x and eq.(A.6) to get: 



i x *(n* 



K' x )Ax 



*x(t) = l T x W x G x (t)W x Pf/(n) x 



(A.7) 



Using the definition of the Laplace transform, f(s) = f^°dtf(t)e st , one finds for the 

moments: 

d n ~ 
(r n )x = (-l) n hml x W' x —G' x (s)W' x Pf/(n) x (A.8) 

There are different ways to obtain expressions for the derivatives needed here. One way 
is to start from the equation for G' x (s), sG' x (s) — E x = W x G' x (s), with E x denoting 
the A^-dimensional unit matrix, to multiply this by powers of s and to differentiate the 
result. This is explained in detail in the Appendices of ref.[15]. This procedure holds for 
arbitrary transition rate matrices, However, in our case, the transition rate matrix W^ 
has a special form because it is the rate matrix for an irreversible system. In particular, 
we can safely assume that there is no eigenvalue A = 0. This is because the matrix II x 
has exactly one vanishing eigenvalue [23]. corresponding to the stationary solution of the 
problem and K.' x consists of at least one non- vanishing matrixelement. Thus, the matrix 
W x has no vanishing eigenvalues and the inverse (W x ) _1 exists. We can thus use the 
form 



G' x (s) = [sE x -W' ; 



G'x(0) 



A' J 



E x + sG' x (0) 



G' x (0) £V \-& x {0) 



(A.9) 



n=0 



with G' x (0) = — (Wjy) l . From this series, one derives 

d n ~ r ~ i "+i 

(-l) n hm— G' x (s) = n\ G x (0) 
s^o as n L 



This yields, using eq.(A.8) 



(r n ) x =n\l x G x (0) Pf/(n) x 



(A.10) 



(A.11) 



Similarly, we o btain fo r the moments of the distribution of consecutive waiting times, again 
using eq.(A.5) and eq.(A.6) along with the two-dimensional Laplace transform f(si,s 2 ) = 



f^dhf^dhfih^e-^e 



'S2<2- 



d 



d n 



r n ) XY = lira \iml Y W' Y —G' Y (s 2 )K x — G' x ( Sl )W' x Pf/(n) x 



2— S-0 Si— 5-0 



ds 



d Sl 



This expression yields with the aid of eq.(A.10|): 
(r n )xY = (n\) 2 l T Y \g' y (Q) K 



x 



Gx(0) Pf/(n) 



x 



(A.12) 



(A.13) 
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From these general equations, eq.(22) and eq.(26) given in the text are derived 



The results for the moments of the turn-over time distribution given in eq.(28) are 
obtained in the following way. One uses the definition (r n ) c = (—1)™ lim s _> 4^V c {s) with 



V c (s) given in eq.(A.3). This expression is rearranged into the form 

V c {s) = l T B W' B G' B (s)K A G' A (s)W' A Pf/(n) A 



by applying the identities (A. 5) and (A. 6). Using eq.( |A.10 ) for the derivatives, after some 
algebra one obtains eq.(28). 



WTDs in the frozen disorder limit 

The limit of vanishing exchange rates , jx;a,a' can be treated analytically because in this 
case one has 

[G' x (t)} a ^ = 5 a , a 'e- kx ^ ; lx , a , a , -> (A.14) 

where kx ]a = ^2 a , kY, a '-,x,a, cf. eq([5]). Using this expression, the WTDs are calculated 
according to eq.(16) and eq.(17) with the result: 



*x(t) = (kx)x 1 E (kx-,*) 2 e- k ^p s l a 

a 

§xAtiM) = (kx)x 1 J2 kY ' ae ~ kYiat2kY ' a ' X ' a,e ~ kx ' a ' tlkx ' 



a'PXa' 



(A.15) 



Furthermore, the moments of $x(t) are given by eq.(34) and eq.(35) and those of 
$x,y(^i,^2) are found to be given by: 



(T n ) X y = (^!)(^)x 1 E^^^'^>^ 



(A.16) 



Eq.(|A.15|) shows that the kinetic scheme is reducible if the rate constants are chosen ac- 

(A.17) 



cording to 

kY,a;X,a' = (kx)^ ' PY,a^Y;a " k X - a < 

and the indices a and a' are independent of each other. 



Appendix B: The phenomenological two-state model 

In this Appendix we give details of the calculations for some expressions relevant for the 
analysis of trajectories. 

Generating functions 



In order to compute F(z,t), one uses eq.(18) and chooses the appropriate expression for 



V, i.e. V c for cycle-counting and V ev in case of event-counting. For a two-state model, the 



corresponding matrices directly follow from eq.(|10|) and eq.(|ll|) and are explicitly given by: 

W' = 



-k A 
k A -k B 



;V C 



k B 
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and 



w 



-k A 
-k B 



j » et; 



fc B 
fc A 



In the expression for the generating function, F(z,t) = l T e^ w +zV )'VP si, /(n), one fur- 
thermore needs the equilibrium populations given in eq.pl). As all matrices are two- 
dimensional, the diagonalization can be performed analytically and the matrix-exponential 
e (w +zV)t can De computed. This way, one finds in case of event-counting: 



" ev\Zj "J 



2w P 



UK 



,(-), 



■,(+), 



:(-)* 



where 



e k ^' l -\- e k ^' t \ + z{l + K) [e k ^' 1 - e fc ™ * 



k f 



Wev = V{1 - K y + AK Z > and fcgO = ^ (1 + K±w. 



Similarly, for cycle-counting one has 

1 



F c (z,t) 



2w c L 



w c e 



K r t 



with 



UK, 



+ e -ki +) A +(l + K) (e- ki ~ h 



kA 



— e 



Kf, t 



y/(l-K) 2 + 4Kz and kf ] = — (1 + K ± w, 



(B.l) 



(B.2) 



Statistics of transitions 

From the generating functions the mean number of transitions and the Mandel parameter 



are calculated as derivatives of the generating function with respect to z as given in eq.(19) 
and one finds 

(k A - k B ^ 2 



W)) 



(ev) 



(1 



-(k A +k B )t 



) + (n)t 



(N(t)) 



(<■) 



2(k A + k B y 

_m_ (e -^» _ 1} + i {n)t 



(B.3) 



where (n) = 2k A kB/(k A + k B \ cf. eq.(31). From these expressions it becomes evident 



that the linear long-time behavior given in eq.(30) is reached after a time on the order 



of (k A + ks)~ l - Also the Q-parameters Q(tY ev ^ and Q(t)^ are time- dependent on the 
same time scale and then approach their constant long-time limits Qoo and Qoo . The 
calculation closely follows the corresponding one of (N(t)), but the resulting expressions 
are not very impressive. The values obtained for the long-time limits are given in the text 



in eq.(32). 



In order to show in which limiting case Poisson statistics is observed, one has to recog- 
nize that this statistics is characterized by an exponentially decaying generating function 
of the form F(z,t) = exp (— (1 — z)nt) with some rate k, yielding {(nt) n /n\)e~ Kt for the 
corresponding probability [T5]. 
It is easy to show that for event-counting one has 



K = l: F ev {t) = e 



-{l-z)k A t _ e -(l-z)k B t 



(B.4) 
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and in case of cycle-counting one finds 

K < 1 : F c (t) = e - {1 - z)kBt and K > 1 : F c (£) = e^ 1 "^* 

Waiting time distributions 



(B.5) 



The WTDs are calculated according to eq.(16) using the fact that here the matrices Kx 



are one- dimensional and consist only of the rates kx and the stationary populations of the 
states are given in eq.Q. As expected, the resulting WTDs are exponential distributions, 
$x{t) — kxe~ kxt , and the corresponding moments are given by (r n )x = n\ky n and in 



v x 



particular (r)x = k x x , cf. eq.(24). Furthermore, ®x,Y{ti,t2) is obtained from eq.(17) and 
one finds that 

$x,y(ii,£ 2 ) = &x(h) • $y(t 2 ) 
which means that this simple kinetic scheme is reducible. 



Appendix C: The TCM 

Definition of the model 

The TCM represents the simplest possible model for dynamic disorder. This model is 
defined as consisting of two configurations a and b in each ensemble A and B, cf. Fig|9j 
This kinetic scheme is irreducible in general. However, there are some limits that give rise 



(A,a) z 



&A. 



a 



o3 






k 



-(B,a) 



B,a 






(A,b)j: 



&A, 






±(B,b) 



&B,b 

Figure 9: Definition of the transition rates in the TCM. We abbreviate the exchange rate 
1x-,b = Jx-,a,b and ^x-,a = Jx-,b,a for X = A, B. Similarly, kx, a and kx,b are abbreviations 
for ky,b-x,a and ky, a x,tn respectively. While the rates kx, a -> ot = a,b depend on the external 
force, k x , a = &x,a(fext)> the exchange rates 7x ;Q are assumed to be independent of f ext . 

to reducible schemes. In the present paper we assume that there is a common equilibrium 
constant, 

ksA^xt) = K(i ext )k Aja (i ext ) a = a, b (C.l) 

meaning that we are essentially treating a fluctuating barrier model[36]. We additionally 
assume that the exchange rates 7x ; b and ^x-,a have the same value, 

lX=lX;a=lX;b ] X = A, B (C.2) 
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In the following, we will not explicitly indicate the force-dependence of the rates and only 
write kx, a and K. Furthermore, it is reasonable to assume that the exchange rates ia 
and 7b show a different dependence on the applied force than the kx,a- For simplicity, we 
assume the that exchange rates do not at all depend on the external force, 7x(f cx t) = 7.x: (0). 
Most quantities of interest are determined by the mean rate and the standard deviation 
for X = A, B: 



k 



X,m 



( k x) = 7T (kx,a + k X ,b) 



1 



kx,d = ox = \/(8k x ) = - {k x , a - k x , 



(C.3) 



Notice that due to eq.(C.l) we have k Bm = KkA, m and k B d — KkA,d- 



Statistics of transitions 

For the TCM one explicitly has the following expression for the relevant matrix of transition 



rates. According to the definition in Fig|9^and eq.(C2) one has: 

/ ~kA,a - 1A 

w = 



\ 



V 



7a k B , a 

lA ~kA,b ~1A k Bjb 

kA,a -k B ,a ~ IB IB 

k A ,b 1b -k B ,b -1b) 



(C.4) 



from which one easily constructs the matrices W and V according the definitions in eq.( 10 ) 



and (11 ). Using the corresponding expressions in the definition of the generating function, 



eq.(18), and the respective derivatives, eq.(19), the mean number of transitions and the 



Q-parameter can be calculated analytically. The long-time limits of the mean number of 
transitions are: 



{N(t)) {ev) = 2(iV(£)) (c) 



K 



KA.mt 



(C.5) 



(1 + K) 

In case of cycle-counting, this result has been given earlier [T9(, [20]. The force-dependence 
of this quantity has been discussed in I. Also the Mandel parameter Q(t) can be calculated 
analytically. The general expressions for the time-dependence of Q(t), however, are not 
very illustrative. Starting from Q(t = 0) = 0, the long-time limit 



Q 



(ev) 

00 



— + Y 

{K + iy x 



K 



K + l 



and Q^ = -2 



K 



X 



K 



(K+l) 2 2 K + l 



(C.6) 



is reached on a time scale that depends strongly on the exchange rate. Here, the deviations 
from Markovian behavior are comprised in the function \i given by 



X 



2k% d [{l + K)\k% r 



k\ d ) + 2k A , m (i A + K-f B 



'l + K)k A , m [{K lA + lB){k 



2 

A,m 



k A,d) + 2 lAlBkA,m] 



(C.7) 



which in the case of a common exchange constant 7 = ja = 1b reduces to the expression 
given in I, x = 2k\ d /(ikA, m )- For small values of the exchange rate, it takes very long 
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until the limit is reached. This can be understood from the limiting behavior for the case 
of static disorder, 'jx — > 0: 



Q(*) (w) = 2Q{t)^ = f A ' d * t (C.8) 



as already found in I [37]. 



Waiting time distributions 

For the WTDs the following bi-exponential expressions are obtained from the definition, 



eq.(16): 



*x{t) = [(k x , m + k 2 X4 )W x cosh (W x t) 

Kx,m,VVx 



+ {lx(k x , m ~ k 2 x , d ) ~ Zk x , m k 2 X4 ) sinh (W x t)] (C.9) 



with W x = yJJx + k 2 x, d - 
Furthermore, one has: 



A<f> AB (tl,t 2 ) = $A,B(tl,i 2 ) - $A(tl)$B(t2 

h 2 (h 2 - h 2 
K A,dK K A,m ^AA 



k 2 , m W A W B 



J^ e -(lA+kA,m)tl e -(jB+kB,m)t2 x 



A,m 

x [(j A + k A , m ) sinh (W A h) - W A cosh (W A ti)] x (CIO) 
x [(7b + k B , m ) sinh (W B t 2 ) - W B cosh (W B t 2 )} 

which is symmetric, A§ A)B (ti,t 2 ) = A$ B ^ A {t 2l ti). 

The moments of <&x{t) are calculated according to eq.(22) and are given by: 



K X,m \ K X,mJ K X,m ~ K X,d + A lxKx,m 

Furthermore, we give the moments of the joint distribution $x,y(^i? £2) i n the form of the 
'deviations from reducibility': 

(Sr n ) X y = (r n )xY ~ (r n )x(r n )Y (C.12) 

For the TCM these moments are symmetric, (5r n )xY — (^ rn )YX- 
With the abbreviation 

Nxy := (k% m - k\ d + 2 lA k A , m ) {K(k% m - k% d ) + 2 TB Au, m ) (C.13) 

one finds: 



kA,d k A k Ad 2 16(k a k Ad )k Ad {^ A + k Am ){^ B + fc s 

h" iv YV and <* r>Afl = — rar ^ 



A,m ^XY ^^A,m ^ XY 

(C.14) 
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